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ABSTRACT 

Numerical solution of DGLAP Q 2 evolution equations is studied for polarized 
parton distributions by using a "brute-force" method. NLO contributions to 
splitting functions are recently calculated, and they are included in our anal- 
ysis. Numerical results in polarized parton distributions and in the structure 
function g\ are shown. In particular, we discuss how numerical accuracy de- 
pends on number of steps in the variable x and in Q 2 . 

Spin-dependent structure functions are measured in polarized lepton-nucleon 
scattering. They depend on two kinematical variables, x and Q 2 . The Q 2 de- 
pendence is calculated within perturbative QCD, and it is described by integrodif- 
ferential equations so called DGLAP equations. Because the Q 2 evolution equations 
are often used in theoretical and experimental studies, it is useful to have a com- 
puter code for solving the equations numerically. Our studies are important, for 
example, in getting optimal parton distributions by analyzing g\ experimental data, 
particularly in obtaining gluon polarization by studying scaling violation of g ± . 

In general, the DGLAP equations are coupled integrodifferential equations: 
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where the variable t is denned by t = — (2//3o) \n[a s (Q 2 )/a s (Ql)], Aqj(x,t) and 
Ag(x, t) are polarized j-flavor quark and gluon distributions, and APy(i) are split- 
ting functions. Each term in Eqs. (la) and (lb) describes the process that a parton 
Pj with the nucleon's momentum fraction y splits into a parton pi with the momen- 
tum fraction x and another parton. The splitting function determines the probability 
of such a splitting process. 

Next-to-leading-order (NLO) splitting functions for spin-dependent parton dis- 
tributions are evaluated recently [1]. Therefore, we can study numerical solution 
of the NLO spin-dependent Q 2 evolution equations [2]. A brute- force method was 
studied for solving spin-independent evolution equations and for those with parton- 
recombination effects [3]. The same method is applied to the spin-dependent case. 
We divide variables into small steps (N x steps in the Bjorken variable In a; and N t 
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Figure 2: N x dependence in xAq s evolution. 

Experimental information on the polarized distributions is, for example, given 
by the structure function g\. In the renormalization scheme MS, the g\ should 
be calculated by the convolution of polarized parton distributions with coefficient 
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We show results of the g\ evolution in 
Fig. 3. The initial distributions are given 
again by the GS set A, and g\ at Q 2 =4 
GeV 2 is calculated by Eqs. (2a) and (2b). 
The #i at Q 2 =200 GeV 2 is calculated first 
by evolving the initial distributions, then by 
taking into account the coefficient functions 
by Eqs. (2a) and (2b). We fix the param- 
eters A^=200 and ^=1000 and show the 

evolved structure functions in the leading Figure 3: Q 2 evolution of xgi . 

order (LO) and in the NLO. Note that the GS distributions, which are obtained 
in the NLO analysis, are also used in the LO input distributions at Q 2 =4 GeV 2 . 
The NLO contributions are conspicuous at small x. Because the gluon polariza- 
tion should be extracted from the scaling violation at small x, our NLO studies are 
important. 

Our numerical analysis are still in progress, particularly in comparison with 
scaling violation data of g\. However, the results indicate that accurate solution is 
obtained in the region 1CT 4 < x < 0.8 by taking more than two-hundred t steps and 
more than one-thousand x steps. 
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